IWOOO&Us 


TDA Progress Report 42-94 


N89 

A ~ 

(ol7V/ 

'1 5p. 


5 / -3,2— 

1018 7 

April-June 1988 



Deep Space Tracking in Local Reference Frames 

R. N. Treuhaft 

Tracking Systems and Applications Section 


A self-calibrating deep space tracking technique is described which can potentially 
produce 2 -nano radian angular spacecraft determinations. The technique uses very long 
baseline interferometric observations of a spacecraft and several radio sources. This article 
first describes the currently employed single-source technique as a parameter estimation 
procedure. Extending the number of parameters and observations leads to the proposed 
local reference frame technique. Station clock , Earth rotation , and tropospheric param- 
eters are estimated along with spacecraft position from the multisource observation 
sequence. The contributions to spacecraft angular uncertainty from system noise , tropo- 
spheric fluctuations , and uncalibrated radio source structure are evaluated. Of these 
experimental errors , radio source structure dominates the determination of the space- 
craft position in the radio reference frame. It is shown, however , that the sensitivity of 
relative spacecraft position accuracies to time-invariant radio source structure effects may 
be on the order of 2 nanoradians. 


I. Introduction 

Deep space angular tracking is routinely performed on 
Voyager with 50- to 100-nanoradian accuracy [1] and will be 
performed on Galileo at the 50-nanoradian level. The analysis 
in this article suggests a tracking strategy for the Deep Space 
Network (DSN) which, in the 30-minute measurement time 
typically allotted for angular tracking, will potentially yield 
better than 2-nanoradian accuracy and will require virtually no 
calibration support external to the signal chains at the stations. 
The self-calibrating technique outlined below capitalizes on 
the extreme sensitivity of the DSN after it is equipped with 
Mark III observing bandwidths and recording hardware [2] . 

The proposed astrometric technique can potentially deter- 
mine angular distances accurate to 10 km near Neptune, or 


2 km near Jupiter. Such a capability would allow the deter- 
mination of a planet’s position in the radio frame by sensing 
gravitational forces long before encounter [3] . l In addition to 
providing navigational benefits, the technique to be described 
would help the DSN realize its potential as a high-precision 
astrometric radio and planetary science tool. This astrometric 
technique would facilitate other tracking-related scientific 
measurement, such as the determination of planetary and 
satellite masses. Additional experiments, such as direct mea- 
surement of relativistic gravitational bending by Jupiter would 
also become possible. 


1 R. N. Treuhaft and J. S. Ulvestad, “Using Gravitational Signatures for 
Target-Relative Angular Tracking During Planetary Approach,” JPL 
Interoffice Memorandum No. 335.3-88-76 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, July 11, 1988. 
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Angular deep space navigation is currently achieved by 
tracking with Very Long Baseline Interferometry (VLBI) 
[4] -[6]. In a technique called Delta Differential One-Way 
Range (ADOR), the geometric delay of a radio source signal 
is differenced from that of a spacecraft signal [1]. Through 
this difference, the angular stability of the radio reference 
frame, determined by the accuracy of the radio source catalog, 
can be transferred to the spacecraft angular measurement. For 
the DSN astrometric radio catalog, measurement errors cur- 
rently limit radio source astrometry to the 10-nanoradian level 
[7] . As measurement techniques and instrumentation improve, 
source structure variability will begin to limit astrometric 
accuracy at the 5-nanoradian level [8] . 

In practice, the current 10-nanoradian stability of the radio 
frame is not completely transferred to the spacecraft angular 
accuracy. Errors in the determination of the spacecraft and 
radio source delays at the time of the ADOR measurement 
limit tracking accuracy to the 20- to 50-nanoradian range, 
which is about 15 to 40 kilometers of projected error at 
Jupiter or 90 to 225 kilometers at Neptune [l]. 2 It has 
been shown that 5- to 10-nanoradian accuracies could be 
achieved with future DSN hardware configurations and im- 
proved calibration support using the ADOR single radio 
source technique. 3 


to which it has been referenced with standard ADOR, it will 
sometimes be in a place where no other close sources are avail- 
able. The local reference frame technique is ideal in situations 
in which a “hole” in the known extragalactic radio sky is 
encountered. 

In the next section, the conventional single radio source 
ADOR strategy is reviewed as a parameter estimation process. 
The local reference frame technique is then developed in Sec- 
tion III as an extension of this process. Section IV contains 
covariance results for the local network technique as com- 
pared to conventional ADOR. Section V summarizes the pro- 
posed technique, outlines possible improvements, and reports 
the status of experimental efforts at validation. 

II. Delta DOR as Parameter Estimation 

In this section, the ADOR process will be formulated as 
parameter estimation from observed VLBI delays. The depen- 
dence of the VLBI delays on spacecraft position and clock 
parameters will be followed by the least squares result for esti- 
mating the parameters and their errors. Extending the depen- 
dence of observed VLBI delays to include station clock rates. 
Earth rotation, and tropospheric effects will then lead to the 
local reference frame technique in the next section. 


The alternative measurement strategy outlined below entails 
interferometric observations of a spacecraft and a number of 
radio sources. This technique of navigation in a local reference 
frame of radio sources can potentially yield 2-nanoradian 
accuracy without clock synchronization, Earth orientation, 
or tropospheric calibration support. The thrust of the local 
reference frame technique is to use the known positions of 
several reference radio sources to determine the magnitude of 
systematic errors affecting the normal ADOR position mea- 
surement. Navigation in a local reference frame is especially 
attractive for outer planet missions in which the spacecraft 
moves by only a few degrees on the sky over a year. For exam- 
ple, Voyager will be in the same 15-square-degree piece of sky 
from late 1987 until Neptune encounter. An additional advan- 
tage of the local reference frame technique is that the radio 
sources do not have to be angularly close to the spacecraft; 
30-degree separations still yield the high accuracy mentioned 
above. When a spacecraft drifts away from the radio source 


The object of angular tracking is to determine the space- 
craft unit vector, which points from the Earth to the space- 
craft. The essential quantities derived from ADOR measure- 
ments are the components of the residual spacecraft unit 
vector. 4 With DSN VLBI, these components are determined 
by observations on nearly orthogonal baselines. The following 
analysis will treat measurements on a single baseline only. In 
that case, the quantity to be determined, A s p , is the compo- 
nent of the residual spacecraft unit vector along the projected 
baseline, which is the baseline projected onto the plane of the 
sky. Given the length of the a priori projected baseline, B p at 
the spacecraft observation epoch, 


In Eq. (1), c is the speed of light, and Ar g is the residual 
spacecraft geometric delay. It is given by 


2 J. B. Thomas, “An Error Analysis for Galileo Angular Position Mea- 
surements With the Block 1 ADOR System ,” J PL Engineering Memo- 
randum No. 335-26 (internal document), Jet Propulsion Laboratory, 
Pasadena, California, November 11, 1981. 

3 R. N. Treuhaft and L. J. Wood, “Revisions in the Differential VLBI 
Error Budget and Applications for Navigation in Future Missions,” 
JPL Interoffice Memorandum No. 335.4-601 (internal document), 

Jet Propulsion Laboratory, Pasadena, California, December 31, 1986. 
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4 The term residual used in the analysis means the difference between 
the actual value of a quantity and its best known a priori value. Resid- 
uals are used both to linearize subsequent least squares procedures and 
to make computation less cumbersome. 
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where $ m is the total baseline vector at the epoch of the space* 
craft measurement, 7 SC is the actual unit vector pointing to the 
spacecraft, and ? msc is the a priori spacecraft unit vector. In 
much of the following analysis, Ar g will be regarded as the 
quantity to be determined from the interferometric observa- 
tions. Geometric delays are related to projected spacecraft 
angles via Eq. (1). The rest of this section is devoted to the 
procedure for estimating Ar g and its error. 


parameter estimate variances due to white observable noise 
are also determined by standard least squares methods. For 
ADOR, with the dependence of observables on parameters 
given by Eqs. (3) and (4), the spacecraft delay and clock epoch 
parameter estimate variances are given by 


Ar g> 


sc,w 


+ a 


rs t w 


0 ) 


The extraction of Ar g from ADOR observations can be 
viewed as a simple example of parameter estimation. In Eqs. 
(3) and (4) below, the observed spacecraft and radio source 
residual delays, A r sc and Ar rj , are expressed in terms of A r g 
and the residual difference between clock epochs at the two 
observing stations, Ar c : 
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where a 2 cvv and o^ sw are the white components of <e 2 c > and 
<e^>, respectively. Note that since white observable noise is 
temporally uncorrelated, it is uncorrelated between the space- 
craft and radio source observations. 

The impact of all correlated errors affecting the delay mea- 
surements can be evaluated in a consider analysis (e.g., [10]). 
This is the formal method for determining how much a given 
error source cancels in the differencing of spacecraft and radio 
source delays prescribed by Eq. (5). The parameter estimate 
variances due to correlated error sources are 


where e sc represents all contributions to observed spacecraft 
delay not included in the first two terms on the right side of 
Eq. (3). Similarly, e rs represents all effects not included in the 
first term on the right of Eq. (4). Since both A r g and A r c are 
delays, they map directly into observed residual delays, and all 
partial derivatives in Eqs. (3) and (4) are unity. Note that a 
residual radio source geometric delay, analogous to Eq. (2), 
does not explicitly appear in Eq. (4). This is because the resid- 
ual delay effects induced by errors in the a priori baseline com- 
ponents and radio source coordinates are included in e sc and 
e rs . The only delay terms which explicitly appear in Eqs. (3) 
and (4) are those which contain parameters to be estimated. 
Generally, both random and systematic effects are included in 
€ sc and e rs , but in the following discussion they are both 
assumed to have zero mean statistics. 


Standard least squares procedures (e.g., [9] ) give the mini- 
mum variance estimates of Ar g and A r c , At and A r c , in 
terms of the observed residual delays: 


At = At -At 

g sc rs 


(5) 


At = At * (6) 

Equation (5) gives the familiar prescription of differencing the 
spacecraft and radio source observed residual delays to esti- 
mate A T g . The estimate of clock epoch given in Eq. (6) is avail- 
able for correcting time tags in other radio metric data. The 
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where oj cc and o} sc are due to correlated components of (ej.> 
and <e£y>, respectively, and p is the correlation coefficient 
between the spacecraft and radio source observations. Familiar 
ADOR rules of thumb, such as “ADOR errors grow propor- 
tionally with spacecraft-radio source separation,” formally 
originate in Eq. (9). Examples of correlated error sources are 
clock rate, Earth rotation, and tropospheric effects. 


III. Local Reference Frame Parameter 
Estimation 

Equations (3) through (10) cast ADOR in the form of a 
parameter estimation procedure. The parameters estimated are 
the residual spacecraft geometric delay and residual differen- 
tial clock epoch. It has been shown that clock rate, Earth 
orientation, and static tropospheric effects are the next domi- 
nant correlated error sources in e sc and e rs (see reference 3). 
(As will be assumed throughout this article, ionospheric errors 
are largely removed by dual-frequency calibration and contrib- 
ute less than 100-picoradian observation errors for the typical 
2300- and 8400- GHz observing frequencies.) The approach 
taken in the reference to reduce ADOR uncertainties from 
the 30- to the 5-nanoradian level was to improve calibration 
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support. Highly accurate earth orientation measurements 
and water vapor radiometry were suggested as the means to 
improving ADOR. An alternative to augmenting calibration 
support is to introduce more parameters into Eqs. (3) and (4) 
and make more radio source observations to solve for them. As 
implied in the introduction, the Mark III instrumentation now 
being installed in the DSN will allow on the order of six highly 
accurate (^30-picosecond) radio source group delay measure- 
ments to be made in the same time currently allotted for one. 


In Eqs. (1 1) and (12) below, Eqs. (3) and (4) are rewritten 
with the extended parameter dependences. The first new 
parameter is r c , which is the delay due to the residual differ- 
ential clock rate between the two stations. The next two 
parameters, 6 and 0. are the magnitudes of rotations about 
two axes perpendicular to the baseline. The rotation axes are 
chosen to be s msc X $ m and Xs mJC Xl m . The next two 
parameters, z 1 and z 2 , are the magnitudes of the static zenith 
tropospheric delay over each of the two stations. The expres- 
sions for all partial derivatives of delay with respect to those 
parameters are given in Appendix A. The local-reference-frame 
extended parameterizations for observed delay are 
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where e^ c and e rs include all observation errors not explicitly 
appearing on the right sides of Eqs. (11) and (12), respectively. 
There are now seven parameters to be estimated. 


It should be noted that the VLBI phase delay rate, which is 
a measure of the short-term rate of change of VLBI delay, is 
extracted along with the delay in the VLBI fringe fitting pro- 
cedure. It is not normally used in standard ADOR because it is 
believed to be dominated by atmospheric fluctuations rather 
than by the geometric signature of the rotating Earth [11] . It 
will be seen, however, that in the local reference frame tech- 
nique, the phase delay rate can be used to constrain singulari- 
ties in the least squares solution for the above seven param- 
eters. The parameter dependences of the observed residual 


spacecraft and radio source delay rates, Ar^ and Af rJ , are 
expressed in Eqs. (13) and (14): 
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The partial derivatives of delay rate with respect to the param- 
eters are given in Appendix A in terms of baseline and source 
coordinates. The new parameter in Eq. (13) is the residual 
spacecraft geometric delay rate, A r g . This parameter, inserted 
for completeness, represents the sum of two contributions to 
the spacecraft’s delay rate: (1) the residual rate induced by 
unmodeled spacecraft motion; and (2) the residual rate in- 
duced by Earth rotation in the presence of a nonzero residual 
spacecraft position vector. In Eqs. (13) and (14), e” c and e rs 
are the contributions to the spacecraft and radio source delay 
rate not explicitly parameterized. 

As will be described in the next section on covariance 
results, the seven parameters in Eq. (1 1) are largely determined 
by delay observations; approximately seven delay observations 
are therefore required in a local reference frame observing 
sequence with the above parameterizations. Only the r g param- 
eter is explicitly determined by the delay rate observation, and 
it is of minimal use in reducing angular position errors. In the 
next section, a data acquisition strategy is evaluated in which 
the spacecraft and five radio sources are observed, with one 
source observed twice. The error in determining residual space- 
craft geometric delay due to system noise, tropospheric fluctu- 
ations, and source structure will be calculated. 


IV. Covariance Results 

A. Observation Configuration for Covariance 
Calculation 

The radio source-spacecraft configuration used to calcu- 
late sample covariance results is shown in Fig. 1. The radio 
sources, which are denoted by Rs in the figure, are from the 
DSN source catalog. Their coordinates in the J2000 reference 
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frame are given in Table 1 [7] . 5 The S in Fig. 1 denotes the 
spacecraft angular position used for the covariance analyses; its 
coordinates are the last entry in Table 1. In these analyses, it 
was assumed that the spacecraft and four of the five radio 
sources were observed only once on the California-Australia 
baseline. One radio source was observed twice to provide the 
seven observations needed to determine the seven parameters 
in Eq. (11) above. The numbers below the letters in Fig, 1 
indicate the temporal sequence of observations. Each observa- 
tion was assumed to last one minute, with two minutes allowed 
for antenna slewing. In Section IVB, the covariance due to 
system noise and tropospheric fluctuations is presented. Sec- 
tion IVC contains the covariance due to source structure 
effects. Although the spacecraft observation was spatially and 
temporally centered in the sample observing sequence, the 
covariance results which follow do not critically depend on 
such a sequence. For example, observing the spacecraft before 
all radio source observations changes the results of the next 
section by only ±1 nanoradian. 

B. Covariance Due to System Noise and 
Tropospheric Fluctuations 

In order to evaluate the spacecraft angular error due to 
white system noise in the local reference frame example of 
Fig. 1, white noise standard deviations of 30 picoseconds were 
assumed for both the radio source and spacecraft observations. 
For the selected radio sources, this assumption is consistent 
with typical Mark III group delay noise values for one-minute 
observations. If it were possible to use the higher-precision 
VLBI phase delays instead of group delays [12] , then typical 
noise standard deviations would be below one picosecond. 
However, this analysis will not assume the availability of the 
phase delay data type. Using standard least squares techniques, 
an equation analogous to Eq. (7) can be derived for the space- 
craft geometric delay error in terms of the system noise con- 
tributions of the seven observations described above. 

The model described in Eqs. (11) through (14) assumes a 
static tropospheric delay over each station. Tropospheric 
fluctuations, predominantly due to the wet component, will 
cause an error in the spacecraft geometric delay. This error has 
been evaluated in a consider analysis using a covariance matrix, 
derived for average DSN conditions, from wet tropospheric 
fluctuation modeling ([11], Eq. [15]). The model in the refer- 
ence assumes that a frozen spatial pattern of turbulence is 
blown across a site by the wind. The model, which is normal- 
ized by WVR data taken near Goldstone, is consistent with a 
daily wet zenith delay fluctuation of 1 cm. 


5 The first source in the table was observed too recently to appear in the 
publication referenced, but its coordinates are derived from the same 
analysis as the other sources in the table. 


Both system noise and tropospheric contributions to space- 
craft angular error are shown in Fig. 2. The ordinate is the 
error in the estimation of A s p of Eq. (1) for the California- 
Australia baseline, which is approximately 10,600 km. This 
angular error is plotted versus hour angle, which is defined 
as the angle between the spacecraft and baseline vectors at 
the middle of the seven-observation sequence. The pattern of 
Fig. 1, therefore, moves across the sky to the west as the hour 
angle decreases. It can be seen that the system noise contribu- 
tion is between 1 and 3 nanoradians. The tropospheric error is 
somewhat smaller. The peak on the right side of the plot is due 
largely to correlations between the first rotation parameter 
and the clock epoch. The location and strength of such peaks 
are completely dependent on the baseline vector and observa- 
tion schedule used. In the example of Fig. 2, if delay rates had 
not been used to constrain the solution, the peak would have 
risen to about 20 nanoradians. The solution could also have 
been controlled by readily achieved a priori constraints on the 
Earth rotation and tropospheric parameters on the order of 
500 nanoradians and 4 cm, respectively. In the covariances in 
this article, the delay rate constraints were chosen to avoid 
placing specific demands on a priori knowledge; but in an 
operational mode, some combination of delay rates and 
a priori constraints could be used to control singularities in the 
determination of geometric delay. 

In order to compare the local reference frame technique to 
a standard ADOR observing sequence, in which Earth rotation 
and static tropospheric effects are not estimated, these effects 
have been evaluated for a spacecraft and single radio source 
observation sequence. The radio source used was the one 
closest to the spacecraft (about 10 degrees away) in Fig. 1, 
P0019+058. Earth rotation uncertainties of 50 nanoradians 
per component and tropospheric zenith delay uncertainties 
of 4 cm were used for the consider analysis. The above Earth 
rotation uncertainty is typical of what will be delivered to the 
DSN for the Galileo and Magellan missions. The tropospheric 
uncertainty of 4 cm reflects the error in wet zenith delay cali- 
bration derived from surface meteorology. Figure 3 shows the 
contributions of these two error sources as a function of hour 
angle. Comparison with the local reference frame system noise 
error of Fig. 2, shown in Fig. 3 as a broken line, shows that 
estimating Earth orientation and tropospheric parameters in 
the multisource observation strategy reduces the tracking error 
due to these effects by an approximate factor of 3 for most 
hour angles. 

C. Covariance Due to Radio Source Structure 

In this subsection, the effect of time-invariant radio source 
structure will be discussed. Time-varying source structure will 
be mentioned in Section V, but no detailed method for mini- 
mizing its effect is given in this article. If the absolute location 
of the spacecraft is required in the radio frame, radio source 
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structure at the 5-nanoradian level imposes a serious limita- 
tion. In the observation configuration of Section A, the abso- 
lute measurement of spacecraft location is degraded to the 
10-nanoradian level due to 5-nanoradian, time-invariant source 
position errors caused by structure effects. Radio source maps 
could be made at regular intervals to calibrate the structure 
effects, but that effort would be costly and in contrast to the 
self-contained nature of the local reference frame technique. 
However, as will be shown below, uncalibrated source struc- 
ture effects can be dramatically reduced for relative angular 
measurements. In many applications, measurement of the 
spacecraft’s change in angular position over some period of 
time can provide valuable navigation information. The stability 
of typical radio source structures over 6-month to 1-year per- 
iods can then be exploited. Because the induced signature of 
constant source structure over these time periods is largely the 
same, it cancels in the relative position measurement. This 
statement is treated rigorously in Appendix B on stationary 
source structure, in which it is shown that maximum error can- 
celing between observations is facilitated by using the same 
observation sequence. 

Using the same observation sequence helps in two ways: 
(1) time-invariant source structure effects on each delay and 
delay rate observation will repeat because the projection of 
the source component vectors onto the baseline will repeat; 
and (2) as shown in Appendix B, observable errors which are 
completely correlated will largely cancel in the differencing of 
sequential spacecraft position measurements if the partial 
derivatives and system noise covariances are similar. The time- 
invariant structure effect will be identically zero if exactly the 
same observation sequence is used. Using similar observation 
sequences means that the spacecraft position cannot change 
substantially or the parameter estimation formalism, which 
leads to the multisource analog of Eq. (5), will propagate the 
structure effect differently from one angular determination to 
the next. This means that, in order to minimize the source 
structure effect and achieve the few-nanoradian relative angu- 
lar accuracy suggested by Fig. 2, the spacecraft cannot move 
significantly in the sky between measurements, which is a 
typical geometry for an outer planet encounter. 

In Fig. 4, the source structure effect on relative angular 
measurements is shown, assuming that the spacecraft moved 
from the S to the indicated points on the two circles between 
measurements. An apparent error of 5 nanoradians in radio 
source position due to structure was assigned for a cos 5 and 
6, where ot is right ascension and 5 is declination, for all radio 
sources. This error is random from source to source, but does 
not vary from epoch to epoch. For the smaller, 3-degree circle, 
which represents approximately 40 days of travel time for 
asymptotic speeds of 10 km/sec, 1- to 3-nanoradian errors in 
spacecraft position estimates are incurred due to source struc- 


ture. For the larger, 6-degree circle, source structure errors 
grow to between 2 and 7 nanoradians. This error growth as the 
spacecraft moves from its initial position is a result of the 
changing least squares equations, as mentioned above. 

V. Summary, Improvements, and 
Demonstrations 

In the absence of radio source structure, the covariance 
analyses presented here show that VLBI observations of 
approximately five radio sources along with a spacecraft can 
yield 2- to 3-nanoradian tracking results. This accuracy, based 
on Mark III radio source measurement noise levels, is attain- 
able through use of the radio source observations to solve for 
clock, Earth orientation, and static tropospheric effects. In 
addition to system noise, the error contributions from tropo- 
spheric fluctuations and radio source structure uncertainties 
have been evaluated in consider analyses. While the tropo- 
spheric fluctuation error contributions are smaller than errors 
due to the assumed VLBI system noise, radio source structure 
can contribute up to 10 nanoradians to the error in locating a 
spacecraft relative to the radio frame. However, this error 
source can be reduced to approximately 2 nanoradians in 
many applications where the relative spacecraft angle over 
small pieces of the sky is of interest. It has been shown else- 
where (see footnote 1) that using a 2-nanoradian differential 
measurement to sense Jupiter’s gravitational field yields target 
relative tracking accuracies competitive with onboard optical 
navigation. However, subnanoradian accuracies will be neces- 
sary to use the same technique at Neptune, for example, and 
improvements to the technique should be considered. 

Improvements to the technique could result from reducing 
the error contributions from system noise, tropospheric fluc- 
tuations, and source structure. As has been mentioned above, 
using the highly accurate phase delay will reduce system noise 
contributions almost two orders of magnitude. This would 
mean that system noise would contribute on the order of 20 
picoradians to the local reference frame technique. Obtaining 
phase delay data is challenging in that it requires resolving 
cycle ambiguities at 8400 GHz or higher radio frequencies. 
Resolving cycle ambiguities is theoretically possible with the 
Mark III system and DSN instrumentation, 6 but this capabil- 
ity has not yet been demonstrated on intercontinental base- 
lines. Even if it is demonstrated, the accuracy inherent in the 
phase delay would be of marginal use if the other error sources 
considered above remained at the quoted levels. 


6 C. D. Edwards, K. M. Liewer, C. S. Jacobs, and R. N. Treuhaft, “An 
Algorithm for Connection of RF Phase in VLBI,” JPL Interoffice 
Memorandum No. 335.3-475 (internal document), Jet Propulsion 
Laboratory, Pasadena, California, November 20, 1986. 


6 


The obvious means of reducing the tropospheric fluctua- 
tion error is by using water vapor radiometers (WVRs) [13] . 
Since typical water vapor fluctuations on 30 -minute time 
scales are on the order of 5 mm [11], the WVR will have to 
track fluctuations at or below the millimeter level to substan- 
tially reduce the fluctuation contribution to the tracking error. 
WVR comparisons with spectral hygrometers, which have been 
conducted in the Advanced Systems Program, suggest the 
possibility of that precision. Note that apparent zenith delay 
biases in the WVR measurement, which have been suggested 
by some existing WVR data [14] , might be tolerated as long as 
they could be parameterized and estimated in place of the 
zenith parameters appearing in Eqs. (11) through (14). 

If WVRs were successfully applied to remove wet tropo- 
spheric fluctuation errors, the dry fluctuations, which are 
estimated to be about 25 percent of the wet on the time scale 
of the observation sequence, would then limit the angular mea- 
surement to approximately 300 picoradians. A possible albeit 
unexplored method of reducing the dry fluctuation effect is 
to measure unmodeled phase fluctuations on Earth-orbiting 
satellites to sense the dry component. Barometric arrays should 
also be considered. Another possibility for improving the angu- 
lar measurement accuracy beyond 300 picoradians is to put an 
interferometric system on an orbiting platform or on the 
moon. In that case, the Earth rotation parameters in the 
system discussed here would be replaced by parameters 
describing the orientation of the platform or the moon. Since 
the variations of these parameters may be very different from 
Earth orientation parameters, the entire local Reference frame 
strategy could conceivably be altered, perhaps ultimately 
reverting back to a single-radio-source ADOR strategy if plat- 
form orientation can be modeled at the 20-picoradian level. 

It has been shown that tracking errors due to time-invariant 
radio source structure can be greatly reduced in relative 
measurements over small pieces of the sky. If tracking rela- 
tive angular position is needed over long periods of time, the 
assumption that the source structure is stationary may not be 
valid. In these cases, it may be possible to parameterize source 
structure evolution 7 and use the accumulated VLBI data to 
estimate and remove the structure effect. The feasibility of 
this technique will be explored as soon as a data base on struc- 
ture effects can be accumulated. 

If, some time in the future, the target ephemeris is known 
to the few-nanoradian level in the radio reference frame, it may 
be desirable to have absolute rather than the differential angu- 


7 J. S. Ulvestad, “Possible Source Structure Effects in IRIS Data,” JPL 
Interoffice Memorandum No. 335.3-88-15 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, February 3, 1988. 


lar measurements discussed above. In other words, determina- 
tions of a spacecraft’s position relative to the radio sources 
may be preferable to determinations of its position relative to 
previous positions. As already mentioned, 10-nanoradian 
absolute errors can be expected from source structure in the 
local reference frame technique. Parameterization of source 
structure may again be a possibility. Another approach for 
decreasing structure-related errors is to use the capability of 
the Mark III system to detect a much larger set of weaker 
radio sources and select only the most pointlike of those 
sources. The trade-off between source structure effects and 
system noise errors will have to be evaluated in this approach. 
It may also be possible to use the structure information 
obtained by observing programs on the Very Long Baseline 
Array [15] to calibrate the structure of sources used in 
local reference frames. 

A potentially important error source which has not been 
considered in this analysis is that of dispersive phase effects 
across the bandpasses of the receiving instrumentation. Because 
the spacecraft signal is a tone and the radio sources are broad- 
band noise, variations in the phase response across the band- 
pass will introduce tracking errors. There may be ways to 
include static bandpass error into the parameter estimation 
scheme. It is also possible that digital front-end electronics 
being explored in the Advanced Systems Program for the DSN 
will allow calibration or elimination of this effect. 

Demonstrations of the local reference frame technique are 
under way to test the assumptions of the covariance analysis 
in this report. Three experiments using the Mark III system 
have been conducted between California and Australia. In the 
first experiment, conducted on March 21, 1988, the gravita- 
tional bend by Jupiter of the ray path of the radio source 
P0201+113 was investigated using the local reference frame 
technique. This gravitational shift in angular position of about 
10 nanoradians simulates a spacecraft motion of about 8 km at 
Jupiter. On April 2 and May 16, when Jupiter was far from the 
ray path, P0201+113 was observed again in the same local 
reference frame of radio sources. About five other groups of 
radio sources were also observed to test the technique. At 
the time of this report, fringes have been found for the first 
two experiments with the Block II correlator at Caltech. How- 
ever, the accuracy of the measurements and the resulting sig- 
nificance of the detection of gravitational bending are as yet 
undetermined. 

If the technique described is validated in a non-real-time 
mode through experiment, an important practical operational 
issue will be the transport of large data volumes in real time. 
As opposed to the current VLBI operational system, which 
relays 500 kbits/ sec from the stations, the Mark III system 
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needed to achieve the derived accuracies will require transpor- 
tation of approximately 50 to 100 Mbits/ sec. Although the 
Advanced Systems demonstrations of the local reference frame 
technique will not be done in real time, a real-time operational 


system would have to employ relay satellites or optical fiber 
links. Data transport rather than analysis will probably be the 
most time-consuming step in producing tracking results for 
navigation. 
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Table 1. Radio source and spacecraft positions used for local 
reference frame covariance analysis 


Radio 

source name 

Right ascension, 
hours/minutes/seconds 

Declination, 

degrees/minutes/seconds 

P0048-09 

0 

50 

41.3180 

-9 

29 

5.216 

P001 9+058 

0 

22 

32.4413 

6 

8 

4.272 

P0106+01 

1 

8 

38.7711 

1 

35 

0.320 

P2320-035 

23 

23 

31.9538 

-3 

17 

5.022 

P2345-16 

23 

48 

2.6085 

-16 

31 

12.019 

Spacecraft 

0 

20 

0.0000 

-5 

0 

0.000 
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Fig. 1. The radio source and spacecraft angular distribution used 
in all covariance analysis results of Section III. In that section, only 
the California-Australia baseline is considered. 




Fig. 3. Standard A DOR covariance results: the standard A DOR 
(single radio source) angular error due to 50-nanoradian Earth 
orientation uncertainty and 4-cm tropospheric delay uncertainty. 
For comparison, the dashed line is the local reference frame system 
noise covariance from Fig. 2. 



HOUR ANGLE, deg 


RIGHT ASCENSION, deg 


Fig. 2. Local reference frame covariance results: the spacecraft 
angular error (As p of Eq. [1 ] ) due to system noise of 30 picoseconds 
per observation and 1-cm daily zenith wet tropospheric fluctua- 
tions, as a function of hour angle. The hour angle is the angle 
between the spacecraft and the baseline at the middle of the 
7-observation sequence. 


Fig. 4. Differential angular position error due to stationary source 
structure: the numbers indicate the error, in nanoradians, in a 
measurement of relative spacecraft position due to time-invariant 
source structure errors. The numbers apply to a differential position 
measurement between the S and the indicated point on the circle. All 
radio sources were assumed to have 5-nanoradian errors in a cos 5 
and in 5 where a is right ascension and 5 is declination. 
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Appendix A 

Partial Derivatives of Delay and Rate 


In the case of standard ADOR, the partial derivatives of 
observed delay with respect to the spacecraft and clock delay 
parameters, and 3 t 5C /3t c , are unity. Below, the 

partial derivatives with respect to all the additional parameters 
in Eq. (11) are expressed in terms of source and baseline com- 
ponents and observation epochs. The partial derivatives of 
delay rate are similarly treated. Since the only differences 
between a spacecraft and radio source observation in this 
treatment are the objects’ angular positions and the observa- 
tion epochs, the partial derivatives below apply to either type 
of observation. The subscripts sc and rs in Eqs. (11) and (12) 
will therefore be ignored. Five partial derivatives, one for each 
parameter added in the local reference frame approach, are 
given below. The derivatives are of delay with respect to 
station-differenced clock rate, two rotation angles, and two 
zenith tropospheric delays. 

The partial derivative of observed delay with respect to 
station-differenced clock rate is 


c 



= f-*n 


(A-l) 


where t is the observation epoch and t Q is the average of all 
spacecraft and radio source observation epochs. 

The first Earth rotation partial derivative in Eq. (11) is as 
follows: 


The partial derivative of observed delay with respect to 
observed station l’s zenith troposphere delay is given by 


3z 1 sin7j 

where y i is the elevation angle at station 1. The error in covari- 
ance results incurred by not using a more realistic zenith par- 
tial derivative is about 0.1 nanoradian. When subnanoradian 
accuracies are attempted, the dry and wet tropospheres, each 
with its own mapping functions [16], [17], will need more 
careful treatment than implied by Eq. (A-3). An expression for 
the partial derivative of geometric delay with respect to station 
2’s zenith troposphere is obtained by substituting 2 for 1 in 
Eq. (A-3) and omitting the minus sign. 

The partial derivatives of delay rate with respect to the 
model parameters of Eqs. (13) and (14) are given below. The 
first three derivatives of rate with respect to spacecraft geo- 
metric rate and clock parameters are 


_3r 

3r 

g 


1 



(A-4) 


— = B , (cos 5 sin a sin 5 - sin 5 cos 6^ sin a ) 

3$ mx\t v c c c' 

+ B , (sin 5 cos 5 cos a - cos 5 cos a sin 6 ) 
m y\t v c c c' 

+ 5 mzlf ( COs6 COs5 c sin ( a c “ a )) ( A ‘ 2 ) 

where, as defined in the text, 0 is the magnitude of a rotation 
about the s msc X axis, where is the a priori baseline 
vector at the spacecraft observation epoch. Note that B mx \ t , 
B mylt , and B mzlt are the a priori baseline components, in a 
space-fixed coordinate system, evaluated at the time t at which 
r is measured. In Eq. (A-2), a, 5, a c , and 5 C are the right ascen- 
sions and declinations of the observed source and the rotation 
axis, respectively. The partial derivative of delay with respect 
to 0, the magnitude of the rotation about B m X t msc X B m , 
is given by an expression identical to Eq. (A-2) with ol c and 
6 replaced by the right ascension and declination of the 

(X, X W Xi L) vect0r - 


where all symbols are defined in the text. The first line of Eq. 
(A-4) applies to the spacecraft rate observation only, as the 
partial of radio source observed rate with respect to spacecraft 
geometric rate is zero. Eq. (A-4) states that delay rate param- 
eters map directly into delay rate observations, and the value 
of the clock epoch parameter does not affect the delay rate 
observations. 

Because an important component of the modeled geometric 
delay rate is due to baseline rotation about the Earth’s spin 
axis, baseline orientation errors cause residual observed delay 
rate. The partial derivatives are found by performing an infini- 
tesimal rotation on the baseline and calculating the difference 
in delay rate between the rotated and original baseline. For a 
rotation about the s m sc X S m axis, with rotation axis right 
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ascension and declination ol c and 8 C , the delay rate partial 
derivative is 

dr o * 
c — = il cos 0 

dd 

X |cosa(i? .cos 5 cos a -B . sin S ) 

L K m z\t c c mx\t c J 

+ sin a (B , cos 5 sin a -B sin 5 Y1 (A-5) 
v m z\t c c my c ; J 


where O is the Earth’s rotation rate in radians/sec. A similar 
expression holds for 3r/30. 

For the partial derivatives of delay rate with respect to 
zenith delay, the troposphere is taken to be static. That is, the 
rate signature is again induced by Earth rotation, which 
causes an elevation angle rate. The partial derivative is found 
numerically by calculating elevation angles at plus and minus 
two seconds from the observation epoch. The rate of change 
of Eq. (A-3) is thereby obtained. 
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Appendix B 

The Effect of Stationary Radio Source Structure on 
Differential Angular Position Determinations 


As noted in Section I VC, time-invariant source structure 
has less effect on differences in angular position than on abso- 
lute angular position in the radio frame. In this appendix, the 
the effect of time-invariant source structure on the angular 
difference determination is evaluated. Figure 4 results from 
the derivation below. The treatment applies equally well to 
radio source position errors, because the source structure 
errors can be represented by unknown shifts in the effective 
radio source positions. For the source structure calculation, it 
will be assumed that the observing schedule for the radio 
sources is identical for the two local reference frame measure- 
ments. The spacecraft will be allowed to move. Its observation 
epoch will be taken to be identical for the two measurements, 
but its position on the sky (and therefore relative to the base- 
line) will be allowed to change. For the following calcula- 
tion, the effect of constant but unknown errors on a general 
parameter estimate, X t , will be considered. This generalized 
notation is consistent with the index manipulations in the 
derivation. For the relative tracking problem in question, X f 
should be thought of as the residual spacecraft geometric delay 
parameter. 


including delay rates in Eqs, (B-l) and (B-2). The F t j and F-j 
coefficients are given by [9] 


F u = ^A T WAY l A T w] i j 
F' u = \{A' T WA'y l A T w] itj 


(B-3) 


The matrix A is defined hy A m ^ n = dT m /dX n . The A matrix 
includes the partial derivatives calculated in the previous appen- 
dix. If the spacecraft has moved between the unprimed and 
primed observations, then the partial derivatives of spacecraft 
delay with respect to model parameters in the A ' matrix will 
be different from those in the A matrix; all other elements of 
A and A ' will be identical to each other. In Eq. (B-3), W is 
the inverse of the modeled observable covariance matrix. It 
has been assumed that the modeled observable covariances, 
such as system noise, do not change between the unprimed 
and primed observations. Although this need not be assumed 
to complete the derivation, the modeled covariances would 
generally be nearly identical for each observation set. 


In the linear least squares formalism, the estimate of the /th 
parameter, X ( , is given by 



y F . A r. 


(B-l) 


where N is the number of observations and the Ar ; s are the 
residual observation delays due solely to source structure 
effects. The F t -s are the least squares coefficients which 
describe the contribution of the / th observation to the / th 
parameter estimate; they are explicitly defined in Eq. (B-3) 
below. Equation (5) in the text is an example of Eq. (B-l) for 
the specific case of ADOR. An estimate of the same param- 
eter, from observations at a later epoch, will be denoted by 


T. F , : .A r i 


(B-2) 


j=l,N 


The A Tj terms are not primed because it is assumed that they 
arise from the same source structure or uncertainty effects for 
both the unprimed and primed observations. Although the 
parameter estimates in this appendix are derived from At- 
delay measurements only, this treatment can be extended by 


In order to calculate the variance of relative X f determina- 
tions due to source structure, the Ar ; ’s will be associated with 
geometric delay errors induced by position uncertainties. 
Below, A Tj is written in terms of parameter residuals A P k : 

(B - 4) 

k k 

where the P k parameters are the right ascensions and declina- 
tions of the observed radio sources. Obviously, the partial 
derivative in Eq. (B-4) vanishes unless P k is the right ascension 
or declination of the radio source being observed in the /th 
VLBI delay measurement. 

The variance of (X t - X ( ) in terms of the covariance of the 
P k source coordinate parameters will now be calculated. Ex- 
pressions for the variance of JQwill be needed. From Eqs. (B-l) 
and (B-4), it is 


var (X.) = 


£ 


9r. 9 t 

. — L M — 
■I bP, Pk,l 9 P, 

k 1 


F. 


i,m 


(B-5) 
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where M p k l is the covariance of the kth and /th source coordi- 
nate errors. Expressions similar to Eq. (B-5) can be derived for 
the variance of X\ and the covariance of X { and X - , leading to 
the final result, shown below in matrix form: 

var (X.-Xj) = (FDM p D T F T ).. + (F'DM p D T F' T ) ii 

-2 (FDM p D T F' T ). . (B-6) 

where D. . = br./bP 
Ui r j 

Equation (B-6) illustrates the parameter error cancellation 
which can arise when parameter estimates are differenced in 


the presence of unknown but repeatable errors. If the observ- 
ing schedules were exactly the same -that is, if the spacecraft 
did not move at all on the sky -then all primed quantities 
would equal all unprimed quantities and Eq. (B-6) would yield 
identically zero. Small changes in spacecraft position make 
small differences between the F and F* and A and^l ' matrices, 
yielding the results of Fig. 4. 8 Note that the crucial step in 
realizing the cancellation was assuming that the Ar ; -’s were 
identical for the unprimed and primed measurements. This is 
equivalent to saying that the source structure or position 
uncertainties did not change from one observation to the next. 


o 

The actual calculation of the results for Fig. 4 employed the gen- 
eralization of Eq. (B-6), which includes both delay and delay rate 
observations. 
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